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1. Introduction 



In a previous study [1] on the vacuum dynamics of non abelian gauge theories we found that 
the deconfinement temperature depends on the strength of an external abehan chromomagnetic 
field. In particular we found that the deconfinement temperature decreases when the strength of 
the applied field is increased and eventually goes to zero. This effect corresponds to the reversible 
Meissner effect in the case of ordinary superconductors therefore we refer to it as "vacuum color 
Meissner effect". This effect could shed light on confinement/deconfinement dynamics, therefore in 
our opinion it is important to test if this effect continues to hold even when switching on fermionic 
degrees of freedom. 

The aim of the present work is to understand if the deconfinement temperature depends on the 
strength of an external abelian chromomagnetic field even in the case of full QCD. 

To this purpose we performed numerical simulations for finite temperature Nf = 2 QCD in an 
external abelian chromomagnetic field. Simulations have been done using APEmille crate in Bari 
and the recently installed computer facilities at INFN apeNEXT Computing Center in Rome. 



2. The method 



To investigate the QCD dynamics in presence of a external background field at finite tempera- 
ture we use the following lattice thermal partition functional [2-4] 



Uk{L,,x)=Ut{0,x)=U^^'(x) 



-(Sw +Sf) 



Uk{L„x)=Ut{0,x)=Ul^\x) 



^Ue'^"^ detM, 



(2.1) 



where S-w is the Wilson action, Sf is the fermion action, and M is the fermionic matrix. The 
spatial links are constrained to values corresponding to the external background field, whereas the 
fermionic fields are not constrained. The relevant quantity is the free energy functional defined as 



-lln. 



(2.2) 



We can evaluate by numerical simulations the derivative of the free energy functional with respect 
to the gauge coupling 



F'(j8) 



dp 



V 



(2.3) 



where the subscripts on the averages indicate the value of the external field. The generic plaquette 
Ufivi^) = U^{x)Uv{x + fi)uj^{x + v)ul{x) contributes to the sum in Eq. (23) if the link Ufx{x) is 



a "dynamical" one, i.e. it is not constrained in the functional integration eq. (|2. Observing that 
/[A*^"'] = at j8 = 0, we may eventually obtain /[A*^"'] from /'[A'^"'] by numerical integration: 



/[A'^'^ViS'. 



(2.4) 
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32 X 8 n =1 
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Figure 1: The derivative of the free energy (blue circles) and the chiral condensate (red circles) vs. j3 on a 
32^ X 8 lattice and bare quark mass anig — 0.075 in correspondence of external field strength next — 1- 

3. Numerical results 



We consider an external abelian chromomagnetic field directed along spatial direction 3 and 
color direction 3. On the lattice this corresponds to constrain the spatial links to 



ur{x)=urix) 



{x) 



exp(/£^) 





exp(- 







)0 

1 



(3.1) 



(3.2) 



where due to toroidal topology of the lattice the field strength gets quantized 

.gH In 
a — = — Wext ,«ext integer. 

In Fig. |l| we display the derivative of the free energy and the chiral condensate at fixed external 
field strength. The peak in the derivative of the free energy correlates to the drop of the chiral 
condensate. 

The deconfinement phase transition is located by looking at the peak of the derivative of the 
free energy. In pure gauge theories we found [1] that the peak position depends on the strength 
of the external chromomagnetic field. Indeed by varying the external field strength we found even 
in the case of full QCD that the phase transition shifts by varying the field strength. In Fig. ^ the 
free energy derivative vs. j8 in correspondence of three values of the field strength is displayed. 
And one can see that the peak of the free energy derivative shifts towards lower values of )8 by 
increasing the field strength. It is also interesting to notice that (see Fig. ||) the value of the chiral 
condensate depends on the strength of the applied chromomagnetic field. 
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Figure 2: The peak position in the derivative of the free energy in correspondence of three values of the 
external field strength. 



4. The critical field strength 

Now we want to repeat the analysis of Ref. [1] in order to ascertain if (i) the deconfinement 
temperature goes to zero as increasing the chromomagnetic field strength (ii) there exists a critical 
field strength beyond which the system is deconfined even at zero temperature. 

The deconfinement temperature can be derived using the values of the critical coupling ob- 
tained by locating the peak of the derivative of the free energy in correspondence of each value of 
the external field strength. 

In our previous study [1] in pure gauge SU(3) we used the string tension as a physical scale 
(see Fig. |). 

Tc 1 

(4.1) 



7^ Lr^MPc)' 



Moreover, using eq. (|3.2|), the field strength is 



y L,-a(/3,) ■ 

The lattice data can be reproduced by a linear fit 



\^Am, ,4.3, 

C7 va va 



with Tc{0) I \fo = 0.643(15). Where our determination for rc(0)/v^ is consistent with the deter- 



minations Tc/\fd = 0.640(15) obtained in the literature without external field [5]. Using eq. (4-.3) 
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Figure 3: The chiral condensate vs. j3 in correspondence of different values of the external field strength. 
The region near the chiral phase transition is displayed in the inset . 



the critical field can now be expressed in units of the string tension 



2.63±0.15. 



(4.4) 



Assuming ^/a = 420 MeV, eq. ( 4^ ) gives for the critical field 

^/^ = (1.104 ±0.063)GeV (4.5) 

corresponding to gHc = 6.26(2) x 10^^ Gauss. 

In the present work we perform a preliminary analysis of the lattice data obtained in full QCD 
with Nf = 2 using the A scale introduced in [6,7] 

1 



A = -fig'){i + c2digy + c4d{g)') 



with 



fig' 



(4.6) 



(4.7) 



fig'- = 1) 

and fig^) the 2-loop jS -function with Nf = 2. 

The data for T^/A versus the external field strength igH)^^^/A are reported in Fig. |5[ In Fig. || 



we display also the pure gauge data analysed using the same scale A given in Eq. ( [4.6[ ). Data 
reported in Fig. (js]) shows that the deconfinement temperature goes to zero, both in quenched and 
unquenched case, in correspondence of a given critical field strength. 

In particular, assuming a Unear dependence on ^/gH as in the quenched case 



TcigH) 
TciO) 



1 



/gH 



(4.8) 
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Figure 4: SU(3) pure gauge. The critical temperature estimated on a 64^ x 8 lattice in units of the string 



tension, eq. (4.1 ), versus the square root of the field strength y/gH in units of the string tension, eq. (4.2) 



Solid line is the linear fit eq. (4.3). In correspondence of zero vertical axis: open circle is Tcj\fo at zero 
external field; full circle is the determination of Td \fo obtained in the literature [5] The green full point is 
the critical field in units of the string tension, eq. (p~4|). 



and using the value Kj ^ = 0.01338 (Eq. 4.3) of Ref. [7]) we found that 

y/^c{Nf = 2) ~ v^^(quenched) ~ l.lGeV. 



(4.9) 



The value for the critical field strength in the present analysis for the quenched case is well con- 
sistent with the value obtained in Ref. [1] using the scale of the string tension and reported here in 



Eq. (^ 



5. Conclusions 

We studied full QCD with Nf = 2 flavors of staggered fermions in presence of a constant chro- 
momagnetic field. We found that, as in the quenched case (see ref. [1]), the critical temperature for 
the deconfinement phase transition depends on the strength of the applied field and eventually goes 
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to zero. The value of the critical field strength is ~ 1 GeV. Moreover numerical results suggest that 
the chiral critical temperature is consistent with the deconfinement temperature and both depend 
on the strength of the external chromomagnetic field. 

Finally our numerical results shows that the value of the chiral condensate increases with the 
strength of the external chromomagnetic field. This is an intriguing effect that deserves further 
investigations. 
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